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Abstract. Using the data of about 33 pb" 1 collected at and around 3.773 GeV with the BES-II detector at 
the BEPC collider, we have studied the exclusive semileptonic decays D + — > K~ir + e + v e , D° -» K 
D + — » K °e + u e and D° — * K*~e + u e . The absolute branching fractions for the decays are measured to 
be BF(D+ -> K~-K + e + u e ) = (3.50 ± 0.75 ± 0.27)%, BF(D° -> K°TY- e + u e ) = (2.61 ± 1.04 ± 0.28)%, 
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1 Introduction 

Semileptonic decays offer access to weak interaction 
matrix elements since the effects of weak and strong inter- 
actions can be separated reasonable well. Measurements 
of branching fractions for exclusive decays of D mesons 
play an important role to develop and to test models of 
their decay mechanisms. 

Earlier theoretical predictions PP0 implied that the 
ratio of the vector to pseudoscalar D meson semileptonic 
decay rates R = r(D ~K* e+v e )/r(D -> ~Ke + v e ) lies in 
the range from 0.9 to 1.2 PJ- E691 Collaboration reported 

a measurement with a lower ratio r{D + — > K e + v e )/ T(D' 
-► K~e + v & ) = 0.45 ± 0.11 0], while MARK-III Collab- 
oration obtained r(D — > K%e + v e )/ T(D — > Ke + v e ) = 
l-O^o 2 El which is consistent with the early theoretical 
predictions. Since then, some controversies evolved con- 
cerning the R ratio from both theoretical predictions and 
experimental measurements. Contrary to the earlier ex- 
pectations, recent theoretical predictions and experimen- 
tal measurements tend to support a smaller value |q][7|- 
A compilation of predictions and experimental values can 
be found in The R ratio is governed by form factors of 
the hadronic currents. So a measurement of R is essential 
to test form factor calculations. 

In this paper, we report direct measurements of the 
branching fractions for the decays D + — > K~n + e + v e 
(throughout this paper, charge conjugation is implied), 

D° -> ~K°iT-e + v ei D+ -> K*°e + v e and D° -> K*~e + u e 
by analyzing the data sample of about 33 pb" 1 8 9 col- 
lected at and around the center-of-mass energy 3.773 GeV 
with the BES-II detector at the BEPC collider. 



2 BES-II detector 

The BES-II is a conventional cylindrical magnetic de- 
tector that is described in detail in Ref. ^Oj- A 12-layer 
Vertex Chamber(VC) surrounding a beryllium beam pipe 
provides input to event trigger, as well as coordinate in- 
formation. A forty-layer main drift chamber (MDC) lo- 
cated just outside the VC yields precise measurements 
of charged particle trajectories with a solid angle cover- 
age of 85% of 47r; it also provides ionization energy loss 
(dE/dx) measurements for particle identification. Momen- 
tum resolution of 1.7%yl + p 2 (p in GeV/c) and dE/dx 
resolution of 8.5% for Bhabha scattering electrons are ob- 
tained for the data taken at y/s = 3.773 GeV. An array of 
48 scintillation counters surrounding the MDC measures 
time of flight (TOF) of charged particles with a resolu- 
tion of about 180 ps for electrons. Outside the TOF coun- 
ters, a 12 radiation length, lead-gas barrel shower counter 
(BSC), operating in limited streamer mode, measures the 
energies of electrons and photons over 80% of the total 
solid angle with an energy resolution of ue/E = 0.22/ \f~E~ 
(E in GeV) and spatial resolutions of = 7.9 mrad and 
<7z = 2.3 cm for electrons. A solenoidal magnet outside 
the BSC provides a 0.4 T magnetic field in the central 
tracking region of the detector. Three double-layer muon 



counters instrument the magnet flux return and serve to 
identify muons with momentum greater than 500 MeV/c. 
They cover 68% of the total solid angle. 



3 Data analysis 

The i/)(3770) resonance is produced in electron-positron 
(e + e~) annihilation at the center-of-mass energy of about 
3.773 GeV. It is believed to decay predominately into 
D°D° and D + D~ pairs. Therefore, if a D meson is fully 
reconstructed (this is called a singly tagged D meson) jHUH], 
the D meson must exist in the system recoiling against the 
singly tagged D meson. In the system recoiling against 
singly tagged D~ and D° mesons, we select semileptonic 
decays D+ -► K-7T + (K*°)e + u e and£>° -> if ti- {K*-)e+v e 
respectively, and measure branching fractions for the de- 
cays directly. 



3.1 Event selection 

In order to ensure the well-measured 3-momentum 
vectors and the reliability of the charged-particle identi- 
fication, all charged tracks are required to be well recon- 
structed in the MDC with good helix fits, and to satisfy 
a geometry cut |cos6*| < 0.85, where 9 is the polar an- 
gle. Each track, except for those from Kg decays, must 
originate from the interaction region, which is defined by 
V xy < 2.0 cm and \V Z \ < 20.0 cm, where V xy and \V Z \ are 
the closest approach of the charged track in the xy-pl&ne 
and z direction. 

Pions and kaons are identified using the dE/dx and 
TOF measurements, with which the combined confidence 
levels (CLtt or CLk) for a pion or kaon hypotheses are 
calculated. A pion candidate is required to have CL n > 
0.001. In order to reduce misidentification, a kaon can- 
didate is required to satisfy CLk > CL^. Electrons are 
identified using the dE/dx, TOF and BSC measurements, 
with which the combined confidence level (CL e ) for the 
electron hypotheses is calculated. An electron candidate 
is required to have CL e > 0.001 and satisfy the relation 
CL e /(CL e + CL K + CLk) > 0.8. 

Neutral kaons are reconstructed through the decay 
Kg — * 7r + 7r~. The difference between the invariant mass 
of 7r + 7r _ combinations and the Kg nominal mass should 
be less than 20 MeV/c 2 . 

Neutral pions are reconstructed through the decay 
7T° — > 77. A good photon candidate must satisfy the fol- 
lowing criteria: (1) the energy deposited in the BSC is 
greater than 70 MeV; (2) the electromagnetic shower starts 
in the first 5 readout layers; (3) the angle between the pho- 
ton and the nearest charged track is greater than 22°; (4) 
the opening angle between the direction of the cluster de- 
velopment and the direction of the photon emission is less 
than 37°. 
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Fig. 1. The distributions of the fitted invariant masses of the 
nKmn (n = 0, l,2;m = 1,2,3,4) combinations in the singly 
tagged D~ modes: (a) K+-R--R-, (b) K -k~ , (c) K°K~, (d) 
K+R-TY-, (e) K°tv-iy-it + , (f) K°iY-ir°, (g) K + -r--r--k°, (h) 
K + ir~Tr~ir~ir + , (i) tt~tt~tt + and (j) sum of the nine modes. 
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Fig. 2. The distributions of the fitted invariant masses of 
the Kmn (m = 1,2,3) combinations in the singly tagged 
D° modes: (a) K + n~, (b) A" + 7r"7r~7r + , (c) K°n + TT~ and (d) 

K+TT-TV . 



3.2 Singly tagged D and D° samples 

The singly tagged D~ and D° samples used in this anal- 
ysis were selected in the previous work 8 9 , we here give 
a brief description for the selection of the singly tagged 
D~ and D samples. 

The singly tagged D~ mesons are reconstructed in 
nine hadronic decay modes of K + tt~tt~, K°it~ , K°K~, 
K+K-JT-, A'°7T-7r-7r+, K°tt-tt q , A'+tt-tt-tt , 
AT + 7r~7r~7r~7r + and 7r~7r~7r + . And the singly tagged D° 
mesons are reconstructed in four hadronic decay modes of 
K + tt~, K + ir~ir~TT + , K°ir + ir~ and K + tt~tt°. 

In order to improve the momentum resolution and the 
ratio of signal to combinatorial background in the invari- 
ant mass spectra, the center-of-mass energy constraint 
kinematic fit is imposed on each of the nKmix (n = 0, 1, 2; 
to = 1,2,3,4) combinations. If there is a K s or 7r° among 
the D daughter particles, an additional constraint kine- 
matic fit will be imposed on the decay K s —>■ ir + ir~ or 
7T° — > 77. Combinations with a kinematic fit probability 
greater than 0.1% are accepted. If more than one combina- 
tion satisfies the criteria in an event, only the combination 
with the largest fit probability is retained. 

The resulting distributions of the fitted invariant masses 
of the nKmn combinations, which are calculated using 
the fitted momentum vectors from the kinematic fit, are 
shown in Fig. 1 and Fig. 2 for the singly tagged D~ and D° 
modes, respectively. A maximum likelihood fit to the mass 
spectrum with a Gaussian function for the D signal and 
a special function jo] |9] to describe the background shape 
yields the observed numbers of the singly tagged D~ and 
D° mesons for each mode. These give the total number of 
the reconstructed singly tagged D mesons, 5321±149±160 
for D~ |E] and 7584 ± 198 ±341 for D° [Hj, where the first 
error is statistical and the second systematic obtained by 
varying the parameterization of the background. 



3.3 Candidates for D + — > K -K + e + v e and 
D° -> K°n-e + p e 

Candidates for the semileptonic decays D + — > K~ ir + e + v e 
and D° — > K iv~e + v e are selected from the surviving 
tracks in the system recoiling against the singly tagged D 
mesons. For the selected candidate events, it is required 
that there should be no extra charged track or isolated 
photon, which has not been used in the reconstruction 
of the singly tagged D mesons. The isolated photon de- 
tected in the BSC should have an energy exceeding 100 
MeV and should satisfy the photon selection criteria de- 
scribed earlier. There are possible hadronic backgrounds 
for each semileptonic decay due to misidentification of 
a charged pion as an electron, for example, the decay 

D + — > AT~7r + 7T + (D° — * K 7r + 7r~) could be misidentified 

as D + — > K~n + e + v e (D° — * K n~e + v e ). However, these 
events can be suppressed by requiring the invariant masses 

of K~Tr + e + (K ir~e + ) combinations to be less than 1.75 
GeV/c 2 . 
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Fig. 3. The distributions of the Umiss for the Monte Carlo 
events of (a) D + — * K~n e v e versus D~ — > K + -k~it~, 
(b) D+ -> K _ 7r + e + z/ e versus D~ -» iT 1 " 7r _ 7r _ 7r , (c) D° -> 

ifVe+^e versus D° -> JsT+tT and (d) D° -» ¥°7r- e + i/ e 
versus D° — » if + 7r~7r°. 
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Fig. 5. The distributions of the fitted invariant masses of the 
nKmn combinations for the events in which the candidates for 
(a) D + — > K~ 7r + e + v e and (b) D° — > K -K~e + v e are observed 
in the system recoiling against the singly tagged D. 
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Fig. 4. The distributions of the Umiss for the selected can- 
didates for (a) D+ -> K~n + e + v e and (b) D° -> K^-K'e+Ve 
from the data. 

In the semileptonic decays, there is one massless neu- 
trino which is undetected. To obtain information about 
missing neutrinos, a kinematic quantity U m i SS = E rn i SS — 
Pmiss is used [HUH], where E m i SS and p m iss are the total 
energy and momentum of all missing particles. To select 
the semileptonic decays, it is required that the candidates 
for the semileptonic decays should have their |i/ m jss,i| < 



3<7[/ mi33 4 , where the <Jjj miss i is the standard deviation of 
the U m iss,i distribution obtained by analyzing the Monte 

Carlo events of D+ -> K-Tt+e+v e (D° -> K°n-e + i> e ) 
versus the ith singly tagged D mode. Fig. 3 shows the dis- 
tributions of the U m i ss i with each peak centered at zero as 
expected, for the Monte Carlo events of D + — + K~ ir + e + v e 
versus D~ — ► K + tt^-k^ , D + — > K~n + e + v e versus D~ — ► 

iC + 7r"7r _ 7r°, D° — ► K n~e + v e versus D° — > K + n~ and 
£>° — > X 7r _ e + ^ e versus 5° — > Jf + 7r~7r°, respectively. 
Fig. 4(a) and Fig. 4(b) respectively show the distribu- 
tions of the U m iss for the selected candidates for D + — > 
/f~7r + e + t/ e and £>° — > if 7r~e + j/ e from the data. 

Figure 5 shows the distributions of the fitted invari- 
ant masses of the nKmir combinations from the events 
in which the candidates for D + — > K ~ir + e + v e and D° — > 
K Ti~e + v e are selected in the system recoiling against the 
nKnm combinations. Fig. 5(a) indicates that an obvious 
signal of D + — > K~ir + e + v e via the D~ tags is observed. 
There are 37 events in the ±3(7m d . mass window of the 
fitted D meson mass , and 18 events in the outside of 
the signal regions, where (Jmd- i s the standard deviation 
of the nKmir distribution for the ith mode. By assuming 
that the distribution of background is flat, 5.6 ±1.4 back- 
ground events in the signal regions are estimated. After 
subtracting the number of background events, 31.4 ± 6.2 
candidates for D + — > K~ir + e + v e are retained. A sim- 
ilar analysis of the events in Fig. 5(b) yields 9.9 ± 3.4 

candidates for D° — > K n~e + v e . There may also be the 
7r + 7r~ combinatorial background. By selecting the events 
in which the invariant masses of the 7r + 7r~ combinations 
on the recoil side of the tags are outside of the Kg mass 
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Fig. 6. The distributions of the invariant masses of (a) K~ir + 
combinations from the selected candidates for D + — > K e + u e 
and (b) K n~ combinations from the selected candidates for 
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Fig. 7. The distribution of the momentum of the electrons 
from the selected candidates for D + — ■» K e + u e , where the 
points with error bars are from the data and histogram is from 
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e v e . 



the Monte Carlo events of D 



K 



window, we estimate that there are 0.6 ± 0.2 background 
events in the candidate events. After subtracting the num- 
ber of background events, 9.3 ± 3.4 candidate events are 
retained. 



3.4 Candidates for D + — » K e + v e and 
D° -> K*-e+v e 

To select the candidates for D + — » K e + u e and 
D° — » K*~e + v e , we calculate the invariant masses of 

K~ir + (K 7r~) combinations from the selected candidates 

for D+ -> K~TT + e + v e (D° -> ~K°ir-e + v e ). Fig. 6(a) and 
Fig. 6(b) show the distributions of the invariant masses of 

K~tt + and K tt~ combinations for the events in which the 
invariant masses of the nKmn combinations are within 
the ±3<7jw d . mass window of the fitted D meson mass 

MD r A clear K signal is observed in Fig. 6(a). Fit- 
ting the K~tt + invariant mass spectrum with a Gaussian 
— *o 

function for the K signal and a S-wave Kit phase space 
background shape, we obtain 29.1 ± 6.6 candidates for 
D + —>■ K e + v e . In the fit, the mass and width of K 
are respectively fixed to 0.8961 GeV/c 2 and 50.7 MeV/c 2 
quoted from PDG |11|. the mass resolution is fixed to 10 
MeV/c 2 determined by Monte Carlo simulation. Using the 
same background shape, a similar fit to the K tt~ invari- 
ant mass spectrum in Fig. 6(b) yields 7.4 ± 3.3 candi- 
dates for D° — > K*~e + v e . After subtracting the number 
of 7r + 7r~ combinatorial background of 0.1 ± 0.1 events, 
7.3 ± 3.3 candidate events are retained. 



However, there are still some K* contaminations from 
other modes of D meson decays or from continuum back- 
ground due to the combinatorial background in the singly 
tagged D signal regions. These K* contaminations must 
be subtracted from the fitted number of the selected can- 
didates for D+ -> K*°e + v e _&nd D° -> K*-e+v e . They 
are estimated by using the D sideband events. The num- 
bers of the events satisfying the selection criteria in the D 
sideband are then normalized to obtain the numbers of the 
background events in the D signal regions. Totally 0.8±0.8 
and 0.7 ± 0.5 background events for D + — > K e + v e and 
D° — > K*~e + v e are obtained respectively. After subtract- 
ing the numbers of the background events, 28.3 ± 6.6 and 

6.6 ± 3.3 candidate events for D + — > K e + i^ e and D° — > 
K*~ e + i> e are respectively retained. 

The distribution of the momentum of the electrons 
from the selected candidates for D + — > K e + v e is shown 
in Fig. 7, where the points with error bars are from the 
data and the histogram is from the Monte Carlo events 
of D + — * K e + v e . The contaminations from hadronic 
and some other semileptonic decays have been subtracted 
from the observed events based on study of the Monte 
Carlo events of e + e~ — > DD. 

3.5 Other backgrounds 

The events from other hadronic or semileptonic decays 
may also satisfy the selection criteria for the semileptonic 
decays and are misidentified as the semileptonic decay 
events. The numbers of the misidentified events have to 
be subtracted from the candidates for the semileptonic de- 
cays. The numbers of the background events are estimated 
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by analyzing the Monte Carlo sample which is about 14 
times larger than the data. The Monte Carlo events are 
generated as e + e~ — + DD, where the D and D mesons are 
set to decay into all possible final states with the branch- 
ing fractions quoted from PDG excluding the decay 
modes under study. The particle trajectories are simulated 
with the GEANT3 based Monte Carlo simulation package 
for the BEST! detector ^2] ■ The number of the events sat- 
isfying the selection criteria is then normalized to the data. 
Monte Carlo study shows that the dominant background 
for D+ -> K-Tr+(K*°)e + is e is from D+ -► lf* /i + ^ 

and the background for D — > K ir~ (K*~)e + v e is from 
D° -► K*-n+. Totally 2.5 ± 0.5, 0.8 ± 0.3, 0.7 ± 0.3 
and 0.2 ± 0.3 background events for D + — > K~n + e + v e , 

D° -> K°TT-e + v e , D+ -> #*°e+i/ e and L>° -> if*-e + i/ e 
are obtained, respectively. After subtracting the numbers 
of the background events, 28.9 ± 6.2, 8.5 ± 3.4, 27.6 ± 6.6 
and 6.4 ± 3.3 signal events for the semileptonic decays 
D+ -> K-TT+e+v e , D" -> Tfn-e+v^ D+ -> K*°e+is e 
and D° — > i4T*~e + f e are obtained. 



where the first error is statistical and the second sys- 
tematic. The systematic error arises mainly from the un- 
certainties in tracking efficiency (—2.0% per track), in 
particle identification (—0.5% per track for charged pion 
or kaon, — 1.0% per track for electron), in photon selec- 
tion (—2.0%), in Kg selection (—1.1%), in U m i SS selec- 
tion (-0.6%), in background subtraction [~(2.5%~9.3%)], 
in Monte Carlo statistics [— (0.8%— 1.4%)], in the num- 
ber of the singly tagged D mesons (—3.0% for D~ and 
—4.5% for D°) and in the fit to the mass spectrum of 

K~ir + or K n~ combination (—1.7% for D + — > K e + i> e 
and —2.5% for D° — > K*~ e + v e ). These uncertainties are 
added in quadrature to obtain the total systematic er- 
ror, yielding -7.6%, -10.8%, -8.0% and -13.7% for the 

semileptonic decays D + — > K~Tr + e + v e , D° — ► K ir~e + v ei 
D + — > K* e + iy e and D° — > if*~e + ^ e , respectively. 

4.3 The ratio of no+^K^e+u.) 

r(D + ->K e+u e ) 



4 Results 

4.1 Monte Carlo efficiency 

The detection efficiencies for the semileptonic decays 



D+ -> K-TT+e+i 



-rrO 



K Tr-e+v e , D+ 



17*° + 



and D° — > X* e + i^ e are estimated to be tD+->K~K+e+i> = 
(15.51±0.12)%,e DO ^o^ e+ ^ = (4.30±0.05)%,e D+ ^.o e . 

= (10.26 ±0.08)% and e D o^ K ,- e +„ e = (2.94 ±0.04)% by 
Monte Carlo simulation, which include the branching frac- 
tions for the intermediate sub-resonance decays. 



4.2 Branching fractions 

The branching fraction for the semileptonic decay 
D —> j (where j = K~7r + e + is e , K ir~e + is ei K e + v e and 



With the measured branching fraction for D + — > 
— *o ■ 

K e ' v e and the previously measured branching fraction 

BF(D+ -> K°e+is e ) = (8.95 ± 1.59 ± 0.67)% by BES Col- 
laboration we obtain the ratio of the vector to pseu- 
doscalar D meson semileptonic decay rates to be 



K* e + v e ) can be determined by 
BF(D - j) 



(1) 



where Njy^j is the number of the signal events for the jth 
mode; is the total number of the singly tagged D 

or D° mesons; er>->j is the detection efficiency for the jth 
mode. Inserting these numbers in Eq. we obtain the 
branching fractions for the semileptonic decays to be 



BF(D + 
BF(D° 



K-TT + e + u e ) = (3.50 ± 0.75 ± 0.27)%, 
K°TT-e + is e ) = (2.61 ± 1.04 ± 0.28)%, 



-*o 



BF(D + — > K e + v e ) = (5.06 ± 1.21 ±0.40)% 



and 



BF(D° 



K*~e + v, 



(2.87 ± 1.48 ±0.39)%, 



r(D+ -> K e+v e ) 
r(D+ ^K°e+v e ) 



0.57±0.17±0.02, 



where the first error is statistical and the second system- 
atic which arises mainly from the uncanceled systematic 
uncertainties including Kg selection (—1.1%), background 
subtraction (—3.5%), Monte Carlo statistics (—1.1%) and 
the fit to the mass spectrum of K~ir + combinations (—1.7%). 



5 Summary 



Using the data of about 33 pb 1 collected around 
3.773 GeV with the BES-II detector at the BEPC collider, 
the absolute branching fractions for the decays D + — > 



- K-jr+e+Ve, D° -> K°TT-e+iy e , D+ 



K*°e+v e and D° 
K ~ix + e + v e ) 



K°ir-e+v e ) = (2.61 ± 
K e^Vp 



K* e + v e are measured to be BF(D + 
(3.50 ± 0.75 ± 0.27)%, BF(D° 

1.04 ± 0.28)%, BF(D+ -> K""e+v e ) = (5.06 ± 1.21 ± 
0.40)% and BF(D° -> K*~e + y e ) = (2.87 ±1.48 ±0.39)%. 

With the measured branching fraction for D + — ► K e + v e 
and the previously measured branching fraction for D + — > 

K e + v e , the ratio of the vector to pseudoscalar D meson 

semileptonic decay rates r(D + — ► K e + v e )/r(D + — ► 

K e + v e ) is determined to be 0.57 ± 0.17 ± 0.02, which is 
in good agreement with theoretical predictions and other 
measurements 00 within error. 
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